Slow relaxation in long-range interacting systems with stochastic dynamics 
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Quasistationary states are long-lived nonequilibrium states, observed in some systems with long- 
range interactions under deterministic Hamiltonian evolution. These intriguing non-Boltzmann 
states relax to equilibrium over times which diverge algebraically with the system size. To test 
the robustness of this phenomenon to nondeterministic dynamical processes, we have generalized 
the paradigmatic model exhibiting such a behavior, the Hamiltonian mean-field model, to include 
energy-conserving stochastic processes. Analysis, based on the Boltzmann equation, a scaling ap- 
proach, and numerical studies, demonstrates that in the long time limit the system relaxes to the 
equilibrium state on time scales which do not diverge algebraically with the system size. Thus, 
quasistationarity takes place only as a crossover phenomenon on times determined by the strength 
of the stochastic process. 
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Systems with long-range interactions have been a sub- 
ject of extensive studies in recent years. In these systems, 
the interparticle potential at large separation, r, decays 
slower than \/r d in d dimensions; for reviews, seellj. Ex- 
amples are self-gravitating systems Q , plasmas [1]^ dipo- 
lar magnets Q , and wave- particle interacting systems Q . 
Long-range interactions lead to nonadditivity, whereby 
thermodynamic quantities scale superlinearly with the 
system size. This often results in unusual features such 
as a negative microcanonical specific heat [6j], inequiva- 
lence of statistical ensembles 0, [|| , and many others [t| . 

Models with long-range interactions often exhibit strik- 
ing dynamical features like slow relaxation [H, [1(3] and 
broken ergodicity [1, [H, 12 1. A very interesting charac- 
teristic feature is the occurrence of long-lived nonequilib- 
rium quasistationary states (QSS) and violent relaxation 



into these states [13|. These states involve a slow relax- 



ation of macroscopic observables over times which diverge 
algebraically with the system size. This suggests that in 
the thermodynamic limit, the system never reaches the 
Boltzmann equilibrium but instead remains trapped in 
the QSS. These intriguing states are usually observed 
under deterministic Hamiltonian evolution. In reality, 
stochastic dynamical moves resulting from coupling ei- 
ther to external environment or to internal degrees of 
freedom are often present. It is thus of interest to ana- 
lyze the robustness of QSS to stochasticity. The aim of 
this work is to explore their existence under stochastic 
dynamics beyond deterministic evolution. 

A prototypical model which allows for detailed studies 
of QSS is the Hamiltonian mean-field (HMF) model. The 
model describes N globally coupled XY spins with the 
Hamiltonian 1141: 
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where &i € [0, 27r) is the phase of the ith spin and pi 
its conjugate momentum. The model represents physi- 



cal systems like gravitational sheet models [15[ and the 
free-electron laser 0. Within a microcanonical ensem- 
ble, the time evolution of the system is described by the 
deterministic Hamilton equations: 
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i y , respectively, are the x and the y 
components of the specific magnetization vector, m — 
^2iLi{cos9i,sin9i)/N . The dynamics, Eq. ([2]), conserves 
energy and momentum. Defining the temperature T as 
twice the specific kinetic energy, the energy per parti- 
cle, e, satisfies the relation e = T/2 + (1 — m 2 )/2, where 
to 2 = to 2 + to 2 and the Boltzmann constant is taken 
to be unity. In equilibrium, the model shows a contin- 
uous transition from a paramagnetic to a ferromagnetic 
phase at the critical energy e c = 3/4 , corresponding to 
the critical temperature T c = 1/2 [3], [l6| . 

In the HMF model, a typical initial state to study re- 
laxation to equilibrium is one which is homogeneous in 
angles and uniform in momenta in some momentum in- 
terval (the "water-bag" initial condition). It is observed 
that, at an energy interval just below e c and under the 
dynamics of Eq. |(5J), the magnetization stays close to its 
initial value over a long time which scales with the sys- 
tem size as N 5 , where 5 > 1 17, Recent studies of 
QSS have revealed features such as anomalous diffusion 
and non-Gaussian velocity distributions [19]. 

Recent studies have invoked a coupling of the HMF 
system to an external heat bath, thereby allowing for 
energy exchange between the two [2(|. These studies 
suggest that QSS may occur, depending on the interplay 
of different time scales underlying the coupling to the 
bath. In particular, it has been shown that the cou- 
pling to the bath results in relaxation to the canoni- 
cal Gibbs-Boltzmann equilibrium state on a time scale 
which does not diverge with the system size. Since qua- 
sistationarity has, so far, been observed only in isolated 
systems with deterministic, energy-conserving dynamics, 
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testing its stability to stochastic energy-conserving pro- 
cesses, where no external bath is involved, would be of 
great interest. 

In this Letter, we address the question of robustness of 
QSS with respect to stochastic dynamics of an isolated 
system within a microcanonical ensemble, where the en- 
ergy is conserved. To this end, we generalize the HMF 
model to include stochastic dynamical moves in addition 
to the deterministic ones, Eq. ([2]). We study our model 
by analyzing the Boltzmann equation for the time evo- 
lution of the phase space density and also by a scaling 
approach and by extensive numerical simulations. We 
provide physical arguments to suggest that the stochastic 
process introduces a cutoff in the relaxation time of the 
QSS. As a result, QSS are maintained only as a crossover 
phenomenon over a characteristic time, which is deter- 
mined by the strength of the stochastic process. In par- 
ticular, at long times, the relaxation time no longer scales 
algebraically with the system size. Consequently, there 
are no quasistationary states. Our scaling form for the 
relaxation time is in very good agreement with results 
from our numerical simulations. 

We now define our generalized HMF model. It fol- 
lows a piecewise deterministic dynamics, whereby the 
Hamiltonian evolution, Eq. ([2]), is randomly interrupted 
by stochastic interparticle collisions that conserve energy 
and momentum. We consider collisions in which only 
the momenta are updated stochastically. Since the mo- 
mentum variable in the HMF model is one-dimensional, 
and there are two conservation laws for momentum 
and energy, one has to resort to three-particle colli- 
sions. Namely, three random particles k) collide and 
their momenta are updated stochastically, (pi,Pj,Pk) — > 
(qi,qj,qk), while conserving energy and momentum and 
keeping the phases unchanged. Thus, the model evolves 
under the following repetitive sequence of events: deter- 
ministic evolution, Eq. ([2]), for a time interval whose 
length is exponentially distributed, followed by a single 
instantaneous sweep of the system for three-particle col- 
lisions, which consists of N 3 collision attempts. 

We proceed by considering the Boltzmann equation of 
our model. In the limit N — > oo, this equation governs 
the time evolution of the single-particle phase space dis- 
tribution f(9,p,t) and is given by 



P- 



df d(v) df 



dl 

dt 



dt 

<9i. _ 

- f(9,p,t)f(8',p',t)f(6",p",t)}, 
R = aS(p +p' +p" -q-q' - q") 



89 dp 

= I dr]R[f(0,q,t)f(e',q',t)f(e",q",t) 



(3) 



(4) 



S[^(p 2 +p' 2 +p" 2 )-l(q 2 + q' 2 + q" 2 y),(5) 

where dr] = dp' dp" dqdq' dq" d9' d9" . In Eq. ©, (v) = 
J dp' J d9' [1 - cos(0 - 6')} f{6',p',t) is the average po- 



tential. Equation §4$ represents the three-body colli- 
sion term, and R is the rate for collisions (p,p',p") — > 
(g, q' , q") that conserve energy and momentum. The con- 
stant a has the dimension of 1/ (time) and sets the scale 
for collisions: On average, there is one collision after ev- 
ery time interval a . The Boltzmann equation with 
similar three-particle collisions in one dimension was con- 
sidered earlier [2l|. We refer to the Boltzmann equation 



with a = as the Vlasov-equation limit Note that 
both the Boltzmann and the Vlasov equations are valid 
for infinite TV and have size-dependent correction terms 
when TV is finite. 

Significant physical insight into the existence of QSS 
can be gained by a direct inspection of Eq. ([3]). In the 
Vlasov limit, corresponding to the deterministic dynam- 
ics, any state which is homogeneous in angles but with 
arbitrary momentum distribution is stationary. In this 
limit, it has been shown that QSS are related to the lin- 
ear stability of the stationary solutions, chosen as the 
initial state 17J, |22j. For example, consider an initial 



state which is homogeneous in angles and uniform in mo- 
menta over [— PO)Po], where po — y/6e — 3. Such a state 
was shown to be linearly stable for energies in the range 
e* = 7/12 < e < e c and unstable for e < e* 
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result, QSS are observed when the energy lies in the range 
e* < e < e c . In a finite system, such states finally relax to 
equilibrium due to finite-size effects which come into play 
over a time ~ N s , where S > 1 )18;]. For example, for 
e = 0.69, numerics gives S ~ 1.7 [17] . At long times, one 
has m(t) — {l/\f~N)e t / NS fort > N 5 , where the prefactor 
accounts for fluctuations in the initial state. For e < e* , 
linear instability results in a faster relaxation towards 
equilibrium as m(t) ~ (l/\> r N)e lt for t ^> I/7. Here, 
7 2 = 6 (7/12 - e) is independent of N [Hj]. Thus, there 
are no QSS for energies below e*. On the other hand, a 
homogeneous state with Gaussian-distributed momenta 
is linearly unstable at any energy below e c , and no QSS 
are observed 0. 

Let us now turn to a discussion of QSS in the gener- 
alized HMF model, i.e., under noisy microcanonical evo- 
lution, in the light of the Boltzmann equation. First, we 
note that, unlike the Vlasov equation, a homogeneous 
state with an arbitrary momentum distribution is not 
stationary under the Boltzmann equation; instead, only 
a Gaussian distribution is stationary. Suppose we start 
with an initial homogeneous state with uniformly dis- 
tributed momenta. Then, under the dynamics, the mo- 
mentum distribution will evolve towards the stationary 
Gaussian distribution. Interestingly, although the mo- 
mentum distribution evolves, the initial 9 distribution 
does not change in time, since for homogeneous 9 distri- 
bution, the p and 9 distributions evolve independently. 

In a finite system, however, there are fluctuations in 
the initial state. These fluctuations make the homoge- 
neous state with Gaussian-distributed momenta linearly 
unstable under the Boltzmann equation at all energies 
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e < e c , as we demonstrate below. This results in a 
fast relaxation towards equilibrium. This equilibration 
mechanism is in contrast to that in short-range systems, 
whereby equilibration is initiated by the formation of 
droplets of critical size followed by coarsening. 

The proof of linear instability of a homogeneous state 
with Gaussian-distributed momenta at ener gies below e c 
was performed in the Vlasov limit in Refs. [17L l22|. We 
have generalized this analysis to the case of the Boltz- 
mann equation 23]. The essential steps of the analysis 
are summarized below, for the simple case of energies just 
below the critical point. 

The stability analysis is carried out by linearizing 
Eq. © about the homogeneous state. We expand 
f(6, P ,t) as f(9,p,t) = /<°>(p)[l + \fW(9,p,t)} with 
= e" p2 / 2T /(27rv / 2^T). Here, since the initial 
angles and momenta are sampled independently accord- 
ing to /'°^(p), fluctuations for finite N make the small 
parameter A of 0{l/\/~N). At long times, the dynam- 
ics is dominated by the eigenmode with the largest 
eigenvalue of the linearized Boltzmann equation, so that 
fW(6,p,t) = f^\p,uj)e^ ke+ut \ Since the average po- 
tential (v) in Eq. ([3]) involves e ±l9 , one needs to consider 
only k = ±1. The coefficients f±± then satisfy [23[ 

2tt 0/(°) 



dp7 (0 V)/i\V^) 



±ipf ^^2if(0) dp 
+ (4tt) 2 J dp'dp"dqdq'dq"Rj^{p')f^{p" 

x [/« (p, u) - f£l (q, U )] = -iu>f£l (p, u). 
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Treating a as a small parameter, we solve the above 
equation perturbatively in a. In the absence of col- 
lisions (a — 0), the above analysis reduces to that 
of the Vlasov equation and to the unperturbed solu- 
tions, namely, the frequencies o/°) and the coefficients 
tt2(q,uW), which are obtained from the analysis in Ref. 
]. In particular, slightly below the critical point e c , the 



unperturbed real frequencies ft^ -* 



,(o) 



are given by 



|Q(°)| sw (2/ % /tF)(T c - T). Thus, in the Vlasov limit, the 
homogeneous state with Gaussian-distributed momenta 
is unstable below the critical energy. To obtain the per- 
turbed frequencies ft to lowest order in a, we now sub- 
stitute the unperturbed solutions into Eq. ([6]). After a 
straightforward but lengthy algebra, one obtains, at an 
energy slightly below the critical point, the perturbed 
frequencies to be given by '\23\ 



ft w |ft (0) |[l + a^], with A 



2tt 3 / 2 
V3 
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This equation suggests that, to leading order in a, the 
frequencies ft are real for energies just below the critical 
value and vanish at the critical point. Thus, a homo- 
geneous state with Gaussian-distributed momenta is lin- 
early unstable under the Boltzmann equation at energies 



just below the critical point and neutrally stable at the 
critical point. 

Armed with the above background, we can now an- 
alyze the evolution of magnetization in our model while 
starting from a water-bag initial condition. The two time 
scales which govern the evolution of the magnetization 
are (i) the scale over which collisions occur, given by a~ 4 , 
and (ii) the scale ~ N s , over which finite-size effects add 
corrections to the Boltzmann equation. The interplay be- 
tween the two time scales may be naturally analyzed by 
invoking a scaling approach, as we demonstrate below. 
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FIG. 1: (Color online) (a) Magnetization vs time for 
iV = 500 at e = 0.69 and for a values (right to 
left) 10" 6 , 10" 5 , 10" 4 , 10" 3 , and 10" 2 . Data averaging has 
been typically over 100 histories. With increasing a, one 
can observe a faster relaxation towards equilibrium. (b) 
ar (a, N)/ In TV vs aN 5 , showing scaling collapse in accor- 
dance with Eq. l|10Jl . Here, e = 0.69. Data averaging varies 
between 5 x 10 4 histories for the smallest system and 100 
histories for the largest one. 

For a -1 < N 5 , and times a" 1 < t < N 5 , the system 
size is effectively infinite and the evolution follows the 
Boltzmann equation. Here, frequent collisions at short 
times drive the momentum distribution towards a Gaus- 
sian. As noted above, until this happens, the initial mag- 
netization does not change in time. Over the time the 
momentum distribution becomes Gaussian, the instabil- 
ity of such a state under the Boltzmann equation leads 
to a fast relaxation towards equilibrium, similar to the 
result for the Vlasov-unstable regime. The asymptotic 
behavior of the magnetization is thus 

m (t)~-L e at ; N s ->f>a-\ (8) 



By requiring that m(t) acquires a value of O(l), the above 
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equation gives the relaxation time rs, determined by the 
stochastic process, as rs ~ InN/a. 

In the opposite limit, a -1 N s , collisions are infre- 
quent, and therefore, the process that drives the momen- 
tum distribution to a Gaussian is delayed. The magne- 
tization stays close to its initial value, and relaxes only 
over the time ~ N s , over which finite-size effects come 
into play. Here, similar to the result for the Vlasov-stable 
regime, the magnetization at late times behaves as 



m(t) 



1 



t/N s 



N 



1 >i> A^ 5 . 



(9) 



This equation gives the relaxation time td, determined 
by the deterministic process, as td ~ N s IniV. 

Interpolating between the above two limits of the time 
scales, one expects the relaxation time r(a, N) to obey 
t- 1 = Tg 1 + t d \ yielding r(a,N) ~ bxN/(a + 1/N 5 ). 
More generally, this suggests a scaling form 



IniV r 
r(a,N) g{aN s ) 



(10) 



where, consistent with Eqs. © and ©, the scaling func- 
tion g{x) behaves as follows: g(x) ~ x for x <C 1, while 
g(x) — > constant for x ^> 1. 

Equation (fT0|) implies that, for fixed N, the relaxation 
time of the water-bag initial state exhibits a crossover, 
from being of order N s In N (corresponding to QSS) for 
a < l/N s to being of order In AT for a > l/N s . This 
brings us to the main conclusion of this work: In the 
presence of collisions, the relaxation at long times does 
not occur over an algebraically growing time scale. This 
implies that, under noisy microcanonical evolution, QSS 
occur only as a crossover phenomenon and are lost in the 
limit of long times. 

To verify the above predictions, in particular, the 
scaling form in Eq. (|10p . we performed extensive nu- 
merical simulations of our model. The Hamilton equa- 
tions, Eq. were integrated by using a symplectic 
fourth-order integrator. In realizing the stochastic pro- 
cess (p,p',p") — > (q,q',q") while conserving the three- 
particle energy E and momentum P, we note that the 
updated momenta lie on a circle formed by the inter- 
section of the plane p + p' + p" = P and the spherical 
surface p 2 + p' 2 + p" 2 = 2E. The radius of this circle 
is given by r = W2E — P 2 /3. The new momenta may 
thus be parametrized in terms of an angle <f> measured 
along this circle, as q = (P/%/3) + r-\/2/3cos0, q' = 
(P/V3) - (r/V6)cos0 - (rlV2)sm<f>,q" = (P/V3) - 
(r / \/6) cos 4>+{r/ \/2) sin <f> [21[ . Stochasticity in updates 
is achieved through choosing the angle <fi uniformly in 
[0,2tt). 

Following the above scheme, typical time evolutions of 
the magnetization in our model for N — 500 and several 
values of a at an energy density e = 0.69 are shown in 
Fig. HJa). The relaxation time r(a,N) is taken as the 



time for the magnetization to reach the fraction 0.8 of 
the final equilibrium value (the result, however, is not 
sensitive to this choice). At e = 0.69, where the equi- 
librium value of the magnetization is ~ 0.3 and <5 ~ 1.7 
17| , we plot ar(a, N) / In A^ vs aN s to check the scaling 
form in Eq. PH|) . Figure [IJb) shows an excellent scal- 
ing collapse over several decades. This is consistent with 
our prediction for QSS as a crossover phenomenon under 
noisy microcanonical dynamics. 

In summary, we have shown that the quasistationarity 
existing with pure deterministic dynamics is lost in the 
presence of stochastic dynamics in the long time limit. It 
occurs only as a crossover phenomenon on times deter- 
mined by the strength of the stochastic process. It would 
be of great interest to explore the general validity of this 
result in other models with long-range interactions. 
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Systems with long-range interactions have been a sub- 
ject of extensive studies in recent years. In these systems, 
the inter-particle potential at large separation, r, decays 
slower than l/r in d dimensions; for reviews, see l|. Ex- 
amples are self-gravitating systems 0, plasmas [||, dipo- 
lar magnets 0] and wave-particle interacting systems Q . 
Long-range interactions lead to non-additivity, whereby 
thermodynamic quantities scale superlinearly with the 
system size. This often results in unusual features such 
as a negative microcanonical specific heat @ , inequiva- 
lence of statistical ensembles 0, Q , and many others @ . 

Models with long-range interactions often exhibit strik- 
ing dynamical features like slow relaxation 0, [l(| and 
broken ergodicity [1, [n], [l2[ . A very interesting charac- 
teristic feature is the occurrence of long-lived nonequi- 
librium quasistationary states (QSS). These states in- 
volve a slow relaxation of macroscopic observables over 
times which diverge algebraically with the system size. 
This suggests that in the thermodynamic limit, the sys- 
tem never reaches the Boltzmann equilibrium, instead 
remains trapped in the QSS. These intriguing states are 
usually observed under deterministic Hamiltonian evo- 
lution. In reality, stochastic dynamical moves resulting 
from coupling either to external environment or to in- 
ternal degrees of freedom are often present. It is thus of 
interest to analyze the robustness of QSS to stochasticity. 
The aim of this work is to explore their existence under 
stochastic dynamics beyond deterministic evolution. 

A prototypical model which allows for detailed studies 
of QSS is the Hamiltonian Mean-Field (HMF) model. 
The model describes N globally coupled XY spins with 
the Hamiltonian IFJ: 



H = 



EPi_ 
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i=l 



1 

2N 



N 
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(6i ~ 0j)} 



(1) 



where 9i € [0, 2ir) is the phase of the i-th spin and pi 
its conjugate momentum. The model represents physi- 
cal systems like gravitational sheet models [l4[ and the 
free-electron laser [E]. Within a microcanonical ensem- 



ble, the time evolution of the system is described by the 
deterministic Hamilton equations: 



dt 



Pi, 



dpi 
dt 



-m r smf 



+ m y cos ( 



(2) 



Here, m x and m y , respectively, are the x and the y 
components of the specific magnetization vector, m — 
J2iLi( C0S ®i> s ' m ®i)/N- The dynamics, Eq. $2$, conserves 
energy and momentum. Defining the temperature T as 
twice the specific kinetic energy, the energy per parti- 
cle, e, satisfies the relation e = T/2 + (1 — m 2 )/2, where 
to 2 = to 2 , + to 2 and the Boltzmann constant is taken 
to be unity. In equilibrium, the model shows a contin- 
uous transition from a paramagnetic to a ferromagnetic 
phase at the critical energy e c = 3/4 , corresponding to 
the critical temperature T c = 1/2 [13] ]. [15j . 

In the HMF model, a typical initial state to study re- 
laxation to equilibrium is one which is homogeneous in 
angles and uniform in momenta in some momentum in- 
terval (the "water-bag" initial condition). It is observed 
that, at an energy interval just below e c and under the 
dynamics of Eq. {5}, the magnetization stays close to its 
initial value over a long time which scales with the system 
size as N s , where S > 1 [TJ, [T3|- Recent studies of QSS 
have revealed many interesting features, e.g., anomalous 
diffusion [l8j], non-Gaussian velocity distributions (l9j . 
and vanishing Lyapunov exponents [19(. 

Recent studies have invoked a coupling of the HMF 
system to an external heat bath, thereby allowing for 
energy exchange between the two [2(j. These studies 
suggest that QSS may occur, depending on the inter- 
play of different timescales underlying the coupling to 
the bath. In particular, it has been shown that the cou- 
pling to the bath results in relaxation to the canonical 
Gibbs-Boltzmann equilibrium state on a timescale which 
does not diverge with the system size. Since quasis- 
tationarity has, so far, been observed only in isolated 
systems with deterministic, energy-conserving dynamics, 
testing its stability to stochastic energy-conserving pro- 
cesses, where no external bath is involved, would be of 
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great interest. 

In this Letter, we address the question of robustness of 
QSS with respect to stochastic dynamics of an isolated 
system within a microcanonical ensemble, where the en- 
ergy is conserved. To this end, we generalize the HMF 
model to include stochastic dynamical moves in addition 
to the deterministic ones, Eq. @. We study our model 
by analyzing the Boltzmann equation for the time evo- 
lution of the phase space density, and also by a scaling 
approach and by extensive numerical simulations. We 
provide physical arguments to suggest that the stochas- 
tic process introduces a cut-off in the relaxation time 
of the QSS. As a result, QSS are maintained only as a 
crossover phenomenon over a characteristic time, which 
is determined by the strength of the stochastic process. 
In particular, at long times, the relaxation time does not 
any more scale algebraically with the system size. Conse- 
quently, there are no quasistationary states. Our scaling 
form for the relaxation time is in very good agreement 
with results from our numerical simulations. 

We now define our generalized HMF model. It follows 
a piecewise deterministic dynamics, whereby the Hamil- 
tonian evolution, Eq. @, is randomly interrupted by 
stochastic inter-particle collisions that conserve energy 
and momentum. We consider collisions in which only 
the momenta are updated stochastically. Since the mo- 
mentum variable in the HMF model is one-dimensional, 
and there are two conservation laws for momentum and 
energy, one has to resort to three-particle collisions. 
Namely, three random particles, (i,j,k), collide and 
their momenta are updated stochastically, (pi,Pj,Pk) — > 
(qi,qj,qk), while conserving energy and momentum and 
keeping the phases unchanged. Thus, the model evolves 
under the following repetitive sequence of events: deter- 
ministic evolution, Eq. ([3]), for a time interval whose 
length is exponentially distributed, followed by a single 
sweep of the system for three-particle collisions, which 
consists of jV 3 collision attempts. 

We proceed by considering the Boltzmann equation of 
our model. In the limit N — > oo, this equation governs 
the time evolution of the single-particle phase space dis- 
tribution f(9,p,t), and is given by: 

'df 
dt 



dj_ df d(v) df 

dt P d9 d6 dp 



(3) 



|Q = j ' d V R[f(e,q,t)f(9',q',t)f(9", q ",t) 



R = aS(p 



- f(0, p ,t)f(e',p',t)f(e",p",t)i 

p'+p"-q-q'-q") 



(4) 



1 



+ p' 2 +p" 2 )-^ + q ,J +q"') ,(5) 



where dn = dp' dp" dqdq' dq" d9' d9" . In Eq. ©, (v) = 
J dp' J dO' [1 - cos(6» - 6')] f(9',p',t) is the average po- 
tential. Equation §4§ represents the three-body colli- 
sion term, and R is the rate for collisions (p,p',p") — > 
(<7, g', q") that conserve energy and momentum. The con- 
stant a has the dimension of 1/time and sets the scale 



for collisions: on average, there is one collision after every 
time interval a -1 . The Boltzmann equation with similar 
three-particle collisions in one dimension was considered 
earlier .21]. We refer to the Boltzmann equation with 
a = as the Vlasov-equation limit [3j. Note that both 
the Boltzmann and the Vlasov equations are valid for in- 
finite N and have size-dependent correction terms when 
N is finite. 

Significant physical insight into the existence of QSS 
can be gained by a direct inspection of Eq. (J3J- In the 
Vlasov limit, corresponding to the deterministic dynam- 
ics, any state which is homogeneous in angles but with 
arbitrary momentum distribution is stationary. In this 
limit, it has been shown that QSS are related to the lin- 
ear stability of the stationary solutions, chosen as the 
initial state [l^, [22J. For example, consider an initial 
state which is homogeneous in angles and uniform in mo- 
menta over [— po,Po]i where po = V6e — 3. Such a state 
was shown to be linearly stable for energies in the range 
e* = 7/12 < e < e c , and unstable for e < e* fll l22j|. 
As a result, QSS are observed when the energy lies in 
the range e* < e < e c . In a finite system, such states fi- 
nally relax to equilibrium due to finite-size effects which 
come into play over a time ~ N s , where S > 1 113. For 
example, for e = 0.69, numerics gives S ~ 1.7 [lfj. AT 
long-times, one has 



m(t) 



1 



J/N" . 



N 



(6) 



where the prefactor accounts for fluctuations in the initial 
state. For e < e* , linear instability results in a faster 
relaxation towards equilibrium as 



1 

t > -, 

7 



(7) 



where 7 2 = 6 (X — e) is independent of N (22j. Thus, 
there are no QSS for energies below e*. On the other 
hand, a homogeneous state with Gaussian-distributed 
momenta is linearly unstable at any energy below e c , and 
no QSS are observed [22| . 

Let us now turn to a discussion of QSS in the gener- 
alized HMF model, i.e., under noisy microcanonical evo- 
lution, in the light of the Boltzmann equation. First, we 
note that, unlike the Vlasov equation, a homogeneous 
state with an arbitrary momentum distribution is not 
stationary under the Boltzmann equation; instead, only 
a Gaussian distribution is stationary. Suppose we start 
with an initial homogeneous state with uniformly dis- 
tributed momenta. Then, under the dynamics, the mo- 
mentum distribution will evolve towards the stationary 
Gaussian distribution. Interestingly, although the mo- 
mentum distribution evolves, the initial ^-distribution 
does not change in time. This is a result of the product- 
measure form of the initial distribution. 

In a finite system, however, there are fluctuations in 
the initial state. These fluctuations make the homoge- 
neous state with Gaussian-distributed momenta linearly 
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unstable under the Boltzmann equation at all energies 
< e c , as we demonstrate below. This results in a fast 
relaxation towards equilibrium. 

The proof of linear instability of a homogeneous state 
with Gaussian-distributed momenta at ener gies below e c 
was performed in the Vlasov limit in Ref . jm |22| . We 
have generalized this analysis to the case of the Boltz- 
mann equation 23]. The essential steps of the analysis 
are summarized below, for the simple case of energies just 
below the critical point. 

The stability analysis is carried out by linearizing 
Eq. ([3]) about the homogeneous state. We expand 
f(6,p,t) as f(9,p,t) = /(°>(p)[l + AfW(8,p,t)} with 
/(o)( p ) = e -p 2 /2T /(27rV2^T). Here, since the initial an- 
gles and momenta are sampled independently according 
to /^(p), fluctuations for finite N make the small pa- 
rameter A of 0(l/v~N). Next, one assumes that, at long 
times, the dynamics is dominated by the eigenmode with 
the largest eigenvalue of the linearized Boltzmann equa- 
tion, so that /W(0,p,i) = 4 1) (p,a;)e 4(fce+ " t ). Since the 
average potential (v) in Eq. Q involves e ±ie , one needs 
to consider only A' = ±1. The coefficients then sat- 
isfy HI 

±»P/i±(P,w) T 

+q(47t) 2 

c(l). 



2tt 



dp'f i0 \p')f ( ±l(p',") 



2i/(°) dp 
dp'dp"dqdq'dq"Rf {Q \p')f^{p") 



x (p, w) - /£V (q, w)] = -iw/i±(p, w). 



(8) 



Treating a as a small parameter, we solve the above 
equation perturbatively in a. In the absence of collisions 
(a = 0), the above analysis reduces to that of the Vlasov 
equation and to the unperturbed solutions, namely, the 
frequencies and the coefficients .f!V(<7, a/ - 1 ), which 
are obtained from the analysis in Ref. 22]. In particular, 
slightly below the critical point e c , the unperturbed real 
frequencies ft (0) = iu/°> are given by « ^L(T C -T). 

Thus, in the Vlasov limit, the homogeneous state with 
Gaussian-distributed momenta is unstable below the crit- 
ical energy. To obtain the perturbed frequencies £1 to 
lowest order in a, we now substitute the unperturbed so- 
lutions into Eq. ([8]) . After a straightforward, but lengthy 
algebra, one obtains, at an energy slightly below the crit- 
ical point, the perturbed frequencies to be given by [23| 



Q « |fi(°)|[l + aA], with A 



2tt 3 / 2 



1 



1 

V5 



(9) 



This equation suggests that, to leading order in a, the 
frequencies fl are real for energies just below the critical 
value, and vanish at the critical point. Thus, a homo- 
geneous state with Gaussian-distributed momenta is lin- 
early unstable under the Boltzmann equation at energies 
just below the critical point, and neutrally stable at the 
critical point. 

Armed with the above background, we can now an- 
alyze the evolution of magnetization in our model while 
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FIG. 1: (Color online) Magnetization vs. time for 
V = 500 at e — 0.69 and for a values (right to left) 
10" 5 , 10" 4 , 10" 3 , 10~ 2 , 10" 1 . Data averaging has been typ- 
ically over hundred histories. With increasing a, one can 
observe a faster relaxation towards equilibrium. 
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FIG. 2: (Color online) ar(a, N)/lnN vs. aN s , showing scal- 
ing collapse in accordance with Eq. (|13p . Here, e = 0.69. 
Data averaging varies between 5 x 10 4 histories for the small- 
est system and 100 histories for the largest one. 



starting from a water-bag initial condition. The two time 
scales which govern the evolution of the magnetization 
are (i) the scale over which collisions occur, given by a -1 , 
and (ii) the scale ~ N s , over which finite-size effects add 
corrections to the Boltzmann equation. The interplay be- 
tween the two timescales may be naturally analyzed by 
invoking a scaling approach, as we demonstrate below. 

For a^ 1 <C N s , and times a" 1 <C t <C N 5 , the system 
size is effectively infinite and the evolution follows the 
Boltzmann equation. Here, frequent collisions at short 
times drive the momentum distribution towards a Gaus- 
sian. As noted above, until this happens, the initial mag- 
netization does not change in time. Over the time the 
momentum distribution becomes Gaussian, the instabil- 
ity of such a state under the Boltzmann equation leads 
to a fast relaxation towards equilibrium, similar to the 
result for the Vlasov-unstable regime given in Eq. ((TJ- 
The asymptotic behavior of the magnetization is thus 



N 6 > t > a 



(10) 



Requiring that m(t) acquires a value of O(l), the above 
equation gives the relaxation time rs, determined by the 
stochastic process, as rs ~ InN/a. 
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In the opposite limit, a -1 3> N s , collisions are infre- 
quent, and therefore, the process that drives the momen- 
tum distribution to a Gaussian is delayed. The magne- 
tization stays close to its initial value, and relaxes only 
over the time ~ N s , over which finite-size effects come 
into play. Here, similar to the result for the Vlasov-stable 
regime in Eq. ([5]), the magnetization at late times be- 
haves as 

m ( t )„ _L e t/N' a - 1 >t>JV J . (11) 

This equation gives the relaxation time td, determined 
by the deterministic process, as td ~ N s InN. 

Interpolating between the above two limits of the 
timescales, one expects the relaxation time r(a, N) to 
obey, r _1 = t^ 1 + , yielding 

hi TV , , 

T(a,N) j— j. 12 

v ' ; a + l/N s K ' 

More generally, Eq. (fl~2"j) suggests a scaling form 

T(a,N)~—g(aN s ), (13) 
a 

where, consistent with Eqs. (fT0|) and (|11[) . the scaling 
function g(x) behaves as follows: g(x) ~ a; for i< 1, 
while g(x) — > constant for i» 1. 

Equation (fT5)) implies that, for fixed AT, the relaxation 
time of the water-bag initial state exhibits a crossover, 
from being of order AT 5 In AT (corresponding to QSS) for 
a < 1/iV 5 , to being of order In TV for a > 1/N S . This 
brings us to the main conclusion of this work: In the 
presence of collisions, the relaxation at long times does 
not occur over an algebraically growing timescale. This 
implies that, under noisy microcanonical evolution, QSS 
occur only as a crossover phenomenon, and are lost in 
the limit of long times. 

To verify the above predictions, in particular, the scal- 
ing form in Eq. (|13[) . we performed extensive numeri- 
cal simulations of our model. The Hamilton equations, 
Eq. were integrated using a symplectic fourth- 

order integrator. In realizing the stochastic process 



(p,p' ,p") — > (q, q', q") while conserving the three-particle 
energy E and momentum P, we note that the updated 
momenta lie on a circle formed by the intersection of 
the plane p + p' + p" = P and the spherical surface 
p2 _|_ p/2 _|_ p//2 _ xhc radius of this circle is given 
by r = W2E — P 2 /3. The new momenta may thus be 
parametrized in terms of an angle <j> measured along this 

circle, as q — + r-J~^ cos <j>, q' = — -^cos</> — 
^= sin g" = -j= - ^cos</> + -^sin</> [21]. Stochas- 
ticity in updates is achieved through choosing the angle 
<j) uniformly in [0, 27r). 

Following the above scheme, typical time evolutions of 
the magnetization in our model for N = 500 and several 
values of a at an energy density e = 0.69 are shown in 
Fig. [TJ The relaxation time r(a, N) is taken as the time 
for the magnetization to reach the fraction 0.8 of the final 
equilibrium value (the result, however, is not sensitive to 
this choice). At e = 0.69, where the equilibrium value 
of the magnetization is ~ 0.3 and 8 ~ 1.7 [Hj], we plot 
ar(a, AT)/ In N vs. aN s to check the scaling form in Eq. 
(fT3l . Figure [2] shows an excellent scaling collapse over 
several decades. This is consistent with our prediction 
for QSS as a crossover phenomenon under noisy micro- 
canonical dynamics. 

In summary, a generalized HMF model, which evolves 
with a piecewise stochastic microcanonical dynamics, is 
introduced and analyzed for the existence of quasistation- 
ary states. It is shown that the quasistationarity existing 
with pure deterministic dynamics is lost in the presence 
of stochastic dynamics in the long time limit. It occurs 
only as a crossover phenomenon on times determined by 
the strength of the stochastic process. It would be of 
great interest to explore the general validity of this re- 
sult in other models with long-range interactions. 
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